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In terms of spinless fermions via the Jordan-Wigner transformation along a snake-like path and 
spin waves modified so as to restore the sublattice symmetry, we investigate static and dynamic 
properties of two-leg antiferromagnetic Heisenberg ladders. The specific heat is finely reproduced 
by the spinless fermions, whereas the magnetic susceptibility is well described by the modified spin 
waves. The nuclear spin-lattice relaxation rate is discussed in detail with particular emphasis on its 
novel field dependence. 
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Spin gaps— the energy gaps in magnetic excitation 
spectra— have been attracting much interest in recent 
years. Haldane [1,2] pioneeringly pointed out possible 
gapped excitations in onc-dimcnsional Heisenberg anti- 
ferromagnets with integral spins. There followed inter- 
esting topics such as spin gaps in a magnetic field [3] and 
mass generation in mixed-spin chains [4,5]. In such cir- 
cumstances, Dagotto et al. [6] stimulated new interest in 
Heisenberg ladder antiferromagnets suggesting another 
electronic mechanism for the gap formation. 

There are a variety of analytic, as well as numeri- 
cal, studies on ladder antiferromagnets. The strong- 
coupling expansion from decoupled dimers [7] and the 
bond-operator representation of rung dimers [8] are the 
very methods to treat two-leg ladders and indeed clarified 
the nature of low- lying excitations. However, they are 
less applicable to thermodynamics. Although resonating 
valence-bond wave functions [9] may be still relevant to 
ladder materials, their variational calculation, even on 
a nearest-neighbor singlet basis [10,11], is practically re- 
stricted to ground-state properties [12,13]. Then we may 
consider making the best use of conventional tools such 
as spin waves and Jordan-Wigner spinless fermions. 

The naivest generalization of the Jordan-Wigner trans- 
formation to coupled spin chains [14] enlightcningly vi- 
sualizes the gap formation with existing interchain in- 
teraction, but the obtained spin gap [15] is less precise 
in the intermediate interchain-coupling regime of our in- 
terest. Recently, performing the Jordan-Wigner trans- 
formation along an elaborately ordered path but simpli- 
fying the mean-field treatment, several authors [16,17] 
obtained a better description of the spin gap for an arbi- 
trary number of legs. We develop their scheme to ther- 
mal and dynamic properties. Another interest in this 
article is to examine a spin-wave picture for spin-gapped 
antiferromagnets [18]. The conventional spin-wave the- 
ory [19,20] applied to Heisenberg single-rung ladders ends 
in diverging sublattice magnetizations and vanishing gap 
[7]. We therefore employ the sublattice-symmetric spin- 
wave theory [21,22], originally applied to square-lattice 
antiferromagnets. Handling Holstcin-Primakoff bosons, 
we make our first attempt to describe ladder antiferro- 



magnets in terms of spin waves. We make interesting 
explorations into ladder systems featuring fermions ver- 
sus bosons in an attempt to provide an intuitive picture 
for understanding the numerical findings and to reveal 
novel quantum phenomena peculiar to low dimensions. 

We study the spin-i antiferromagnetic Heisenberg 
model on two-leg ladders: 
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which has isotropic nearest-neighbor exchange interac- 
tions along the chains (J) and along the rungs (J'). In 
the two extreme cases of J' = and of J = 0, the 
low-energy properties are well known. The decoupled 
chains are critical, whereas the decoupled rung dimers 
are massive. The intermediate region with J' ~ J 
is thus interesting, but most of model materials lie in 
the regime of large J', such as (^(CsH^^^CLt with 
J' ~ 5J [23], Cu 2 (C5Hi2N 2 )2Br4 with J' ~ 7J [24], and 
(C 5 Hi 2 N) 2 CuBr4 with J' ~ 3.5 J [25]. Hence much at- 
tention is paid to a copper oxide SrCu2 03, comprising 
two- leg ladders with J' ~ J [26]. We consider the case 
of J 1 = J unless otherwise noted. 

It is along a snake-like path [16,17], = (1,1) — > 
(2,1) -> (2,2) -► (1,2) -> (1,3) ->■■••, that we align 
spinless fermions (SFs). When we introduce renum- 
bered spin operators Si.j — Si j (Sjj) for odd (even) 



j's, where i = 3 



the SFs are created as cj , = 



The fermionic Hamiltonian is given by 
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and we search for its mean-field solutions. Setting the 
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equal to 1/2 under zero magne- 



tization and defining a unitary transformation 
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we can diagonalize the Hamiltonian as ^ fe |efc|(a| k ai,k — 
a t fe a 2,fc), where 



Compared with the conventional spin-wave theory, where 
spins on one sublatticc point predominantly up, while 
those on the other predominantly down, the modified 
spin waves (MSWs) restore the sublattice symmetry. If 
we define the Bogoliubov transformation as 
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with xo = (4jCij), Xi = (4,jCi,j + i), and X 2 = 

(c\ jC2j+i) to be self-consistently determined at each 
temperature in an approximation of the Hartree-Fock 
type, while within the Hartree-level approximation, e k — 
J 1 '/2 + iJsin/c, suggesting a spin gap J'/2. 

A bosonic approach starts from the Holstein-Primakoff 
transformation: = (2S — a\ jai.jY^ai.j , Sfj = S — 
a\ jO,i,j for even i+j's, while S^. = a\^{2S—a\-a2,jY^ 2 , 
Sfj = —S+a\ jCL2,j for odd i+j's. Then the Hamiltonian 

is expanded as H — Y^hL-2 wnere is the 0{1/S l ) 
term and is taken into consideration up to i = 0. H-i 
and 7io contain bilinear and biquadratic terms with re- 
spect to the bosonic operators and describe the linear 
and interacting spin waves, respectively [27]. In order to 
avoid quantum as well as thermal divergence of the sub- 
lattice magnetizations, the Bogoliubov transformation is 
carried out subject to the constraint that the sublattice 
magnetizations be zero [21,22]: 
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FIG. 1. The ground-state energy (a) and the spin gap (b) 
as functions of J' /J, and the dispersion relation of spin-triplet 
excitations (c). 
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where A is the Lagrange multiplier due to the constraint 
(5) and e k is the MSW dispersion relation: 

e fe = y/[(2J + J')S + 2JX} 2 - [(J' -2J cos k)S} 2 . (8) 

The optimum thermal distribution functions of 
the MSWs are given by (f3j k f3i,k) = n»,k = 
[ e (£ fc - 9m f/co SW r)ABT _ X ]-i_ We consi( ier H as a per- 
turbation to H-i [29]. 

First we calculate the ground-state energy E g and the 
spin gap A by the linear MSWs (LMSWs), the interact- 
ing MSWs (IMSWs), the Hartree-level SFs (HSFs), and 
the Hartree-Fock-level SFs (HFSFs), and in Fig. 1 show 
them all together with numerical (Exact) findings [7] and 
strong-coupling-expansion (SCE) results [30], E g /NJ = 
-3r/4-3/8r + 0(r- 2 ) and A/ J = r-l + l/2r + 0(r- 2 ), 
where r = J' /J. Unless J' is sufficiently large, the 
ground-state energy is better described by the MSWs, 
while the spin gap by the SFs. When we focus on the in- 
teresting point J' = J, the IMSWs give the best estimate 
of the energy as E g /NJ = —1.144, compared with the 
numerically exact value —1.156, whereas the HSFs give 
that of the gap as A/ J = 1/2, compared with 0.50(1). 
The HSF description is no more precise once the system 
moves away from this point, but the HFSF description re- 
mains quantitative in a wider range of J'/ J. The MSWs 
are trivially less relevant to decoupled clusters, where 
no "wave" can spread over the system. The dispersion 
relation is also shown in Fig. 1. The band bottom is 
well reproduced by the SFs, while the band width by the 
MSWs. Since the MSWs considerably underestimate the 
spin gap, any thermally activated behavior at low tem- 
peratures should be observed through the SFs. 

Secondly we calculate the thermal properties and com- 
pare them with quantum transfer-matrix (QTM) calcu- 
lations [31]. The specific heat C is shown in Fig. 2(a). 
The SFs well reproduce the Schottky peak, while the 
MSWs lose their validity with increasing temperature. 
The HSFs (HFSFs) underestimate the high-temperature 
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behavior C/Nk B = (3k B T/4J) 2 + 0[(k B T/ J) 4 ] by a fac- 
tor 2/3 (4/9) but complement numerical tools very well 
at low temperatures. If we approximate the dispersion 
relation of the low-lying excitations as 



E k = A + Ja(k - tt) 2 , 



(9) 



where A = J'/2 and a = J / J' in the Hartree approxima- 
tion, while A = |(l/2-xo)J'-(l/2-Xi-2x?-X2)J| and 
a = (l/2- X o)[(l/2-Xi-2Xi-X2)J' + 4(x2 + 2xi)J]/2Zi 

in that of the Hartree-Fock type, the SFs illuminate the 
low-temperature behavior as 
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Considering the difficulty of grand canonical sampling at 
low temperatures by numerical tools such as quantum 
Monte Carlo and density-matrix renormalization group, 
the SFs can play an effective role in our explorations. Fig- 
ure 2(b) shows the magnetic susceptibility \- The SFs 
overestimate the round peak at intermediate tempera- 
tures but again work well at low temperatures revealing 
the initial increase as 
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The MSWs are good at reproducing the overall behav- 
ior. All the calculations converge into the paramagnetic 
susceptibility 2S(S + l)/3k B T at high temperatures. 

Lastly we investigate the nuclear spin-lattice relax- 
ation rate in terms of the SFs, which are much better 
than the MSWs at describing the spin gap and there- 
fore the low-temperature properties. Considering the 
electronic-nuclear energy-conservation requirement, the 
Raman process plays a leading role in the relaxation, 
which is formulated as 
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where A is the hyperfine coupling constant between the 
nuclear and electronic spins, wn = JnH is the Larmor 
frequency of the nuclei with 7n being the gyromagnetic 
ratio, and the summation ^ n is taken over all the elec- 
tronic eigenstates \n) with energy E n . Taking account of 
the significant difference between the electronic and nu- 
clear energy scales (Tiwn ^ W~ 5 J) and assuming a rea- 
sonable temperature range k B T < A for SrCu203 with 
A/k B ~ 420 K [26], the relaxation rate is represented as 
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FIG. 2. The specific heat (a) and the magnetic susceptibil- 



ity (b) as functions of temperature. 
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FIG. 3. Temperature (a) and field (b) dependences of the 
nuclear spin-lattice relaxation rate calculated by the HFSFs. 
Numerical estimates of Eq. (13) (solid lines) are compared 
with the analytic expression (14) (dotted lines). 
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where we have again employed the approximate disper- 
sion (9). The term 71^(1 — n^k) is the consequence of the 
principle of detailed balancing in a fermionic ensemble. 
At moderate fields and temperatures, k B T <C A — g[i B H, 
Eq. (13) can be further calculated analytically as 



where K is the modified Bessel function of the second 
kind and behaves as Kq(x) ~ 0.11593— lna; for < x « 1. 
Equation (13), together with its approximate expression 
(14), is plotted in Fig. 3. 

At low temperatures, 1/Ti also exhibits an increase 
of the activation type but with logarithmic correction, 
which is much weaker than the power correction in the 
case of the susceptibility. Such a pure spin-gap-activated 
temperature dependence of 1/Ti, which was pointed out 
in a sophisticated numerical work [31] as well, may be ob- 
served at sufficiently low temperatures but is not yet ver- 
ified for ladder antiferromagnets such as SrCu203 owing 
to the magnetic impurities masking the intrinsic proper- 
ties. Equation (13) deviates from the simple spin-gap- 
activated behavior (14) for k B T > 0.2 J ~ 0.4Z\, which is 
interestingly consistent with the criterion for the break- 
down of the simple two-magnon scheme observed in a 
spin-1 Haldane-gap antiferromagnet AgVP2S6 [32]. In 
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such a high-temperature range, we have to fully incorpo- 
rate particle collisions into the calculation in order to 
describe the crossover to diffusive behavior, where an 
enhanced activation gap may be observed due to the 
temperature-dependent diffusion constant [33]. In the 
low-temperature range /cbT < 0.2J, on the other hand, 
the relaxation rate should strictly follow the expression 
(14): With increasing field, 1/Ti first decreases logarith- 
mically and then increases exponentially. The initial log- 
arithmic behavior comes from the Van Hove singularity 
peculiar to one-dimensional energy spectra and may arise 
from a nonlinear dispersion relation at the band bottom 
in more general. Therefore, besides spin-gapped antifer- 
romagnets, one-dimensional ferro- and fcrrimagnets may 
exhibit similar field dependence [34]. The present log- 
arithmic field dependence at low temperatures should 
and could be distinguished from the 1/VH or ln(l/if) 
dependence of diffusion-dominated dynamics [35,36] at 
high temperatures. The following exponential increase 
originates in the spin-1 excited state lowering in en- 
ergy with increasing field. If we consider the strong- 
correlation compound SrCu 2 3 with J/fc B — 840 K [26], 
the minimum is supposed to appear at H ~ 20 T. In 
the case of molecular-based ladder systems [23-25] with 
much smaller exchange interactions, 1/Ti should reach a 
minimum at much smaller fields. We expect that such 
intrinsic features of the low-energy dynamics of ladder 
antifcrromagnets will be found experimentally. 

We have demonstrated new schemes of investigating 
Heisenberg two-leg ladder antiferromagnets. Although 
both fermionic and bosonic approaches possess advan- 
tages of their own, we stress the SFs along the snake-like 
path as one of the best languages for the present system. 
One of the fatal weak points in the MSW description 
is nonvanishing specific heat at high temperatures. The 
endlessly increasing energy with increasing temperature 
is because of the temperature-dependent energy spec- 
trum (8), where A, playing the role of the chemical po- 
tential, turns out a monotonically increasing function of 
temperature. Such a difficulty can be overcome in ferri- 
magnets exhibiting noncompensating sublattice magneti- 
zations [37,38] but plagues bosonic approaches applied to 
antiferromagnets. The SFs well reproduce the spin gap 
and thus reliably describe the low-temperature proper- 
ties. One of the most interesting findings is the novel 
field dependence of\/T\. Relevant relaxation-time mea- 
surements are strongly encouraged. 
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